THE GAUSSIAN FREE FIELD AND HADAMARD'S VARIATIONAL FORMULA 
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Abstract. We relate the Gaussian free field on a planar domain to the variational formula of 
Hadamard which explains the change of the Green function under a perturbation of the domain. 
This is accomplished by means of a natural integral operator - called the Hadamard operator - 
associated with a given flow of growing domains. 



1. Introduction 

The Gaussian free field (GFFq) on a planar domain O is an analogue of the classical Brownian 
motion with two temporal dimensions, which frequently appears as the continuum limit of 
discrete random processes and plays an important role in statistical physics (see [11 1 for a 
general survey). It can be defined, at least formally, as a random linear combination of an 
orthonormal basis ipj (/ = 1, 2, 3, . . .) of the Sobolev space Hj(0), 

where £.j {j - 1,2,3, . . .) is a sequence of independent standard real-valued Gaussian variables 
(i.e., taken from N(0, 1)). The space Hg(Q) consists of real-valued functions and is equipped 
with the usual Dirichlet inner product. The resulting stochastic field ^ is conf ormally invariant 
with covariance structure determined by the Green function of the domain; thus the GFFq also 
exhibits interesting connections to complex analysis and potential theory. 

The purpose of this note is to relate the Gaussian free field to the classical variational formula 
due to Hadamard |5| which describes the change of the Green function under a perturbation 
of the boundary of the imderlying domain. We interpret Hadamard's formula in terms of a 
natural and concrete integral operator - called by us the Hadamard operator - associated with 
a given flow of smooth domains. It supplies an isometric isomorphism L^(0) Hj(Q), and, 
consequently, it can be used to produce the corresponding GFFq field from the standard white 
noise field. This construction has certain appealing features since it enables us to view GFFq 
as formed by adding independent infinitesimal harmonic fields whose boundary amplitudes 
correspond to suitably weighted white noise processes. 

The paper is organized as follows. In Section|2l we recall the definition of the space (Q) and 
some related prerequisites. In Section|3) we first describe a version of Hadamard's variational 
formula, and then proceed to introduce the associated Hadamard operator and analyse its 
properties. Section|4]contains a brief account on the white noise and Gaussian free fields from 
the point of view of Gaussian Hilbert spaces. Finally, in Section |5) we use the Hadamard 
operator to construct the Gaussian free field and derive some natural consequences. 

Date: 8 February 2012. 

2010 Mathematics Subject Classification. 60K35, 30C70. 

Both authors were supported by the Goran Gustafsson Formdation (KVA) and Vetenskapsradet (VR). The second 
author was in addition supported by the Academy of Finland, project 136785. 

1 



2 



HAAKAN HEDENMALM AND PEKKA J. NIEMINEN 



2. Preliminaries 

In this section, Q is a finitely connected bounded domain in the complex plane C with a 
C^-smooth boundary. We use the notation 

A - ^ 



dy^' ' \dx' dy 

for the Laplacian, the gradient (nabla) operator, and the area element, respectively; here, 
z = X + iy is the decomposition of z e C into real and imaginary parts. 

2.1. The Green function and Poisson kernel. The Dirichlet Green fimction Gq, in Q solves, 
for w e Q, 

{-^zGa{z, w) = 6a,(z), z e Q, 
\Gq{z,w)-Q, z e ^Q, 

where A2 indicates that the Laplacian is taken with respect to z. It has the symmetry property 
Ga{z,w) = Gn{w,z). We write Ga for the associated integral operator 



Go/(z):= r Ga(z,w)f{w)dA{w), 
Jn 



z e Q, 



extended to vanish in C \ Q. Then, in the sense of distributions, —AGaf = / in Q, so we may 
write Gq = [-A]-i. 

The Poisson kernel Pq is obtained from 

Po(z, Q := ^^Mz, C), z e Q, C e da, 

where the unit normal derivative is taken in the inward direction. For C £ dCl, we extend 
Pq(z, C) to (area-a.e.) z e C by setting it equal to in C \ cl[Q], where cl[Q] denotes the closure 
of Q. The associated integral operator is 

Pa/(z):= r PQ(z,Q/(Qds(0, z e Q, 
Jda 

where ds refers to the arc-length measure. It furnishes the harmonic extension of / defined on 
to the interior Q. 

2.2. The space Hj(Q). The Sobolev space Hj(Q), alternatively denoted by Wg'^(Q), can be 
defined as the Hilbert space completion of (real- valued) C^(Q) under the Dirichlet inner 
product 



</,g>Q,v- V/-VgdA, 
Jo 

where the dot is the inner product in IR^. The norm of the space is given by ||/||q,v (///)q v' 
As all functions f,g ^ Hg{Q) vanish on the boundary, we may write, using integration by parts, 

(2.1) if, g}ay = (-A/, = {[-A]''^f, [-An)n, 

where (-, •)q is the usual inner product in L^(Q). It follows that the operator [— A]^^^ supplies 
an isometric isomorphism from Hj(Q) onto L^(Q), and its inverse is given by [— A]~^^^. As 
is standard, the powers of —A are understood in terms of eigenfunction expansions with 
vanishing Dirichlet boimdary data (cf., e.g., |^ Sec. 6.5.1]). 
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3. HaDAMARd's FORMULA AND AN ASSOCIATED OPERATOR 



3.1. Hadamard's variational formula. We consider a family of bounded planar domains Q(f) 
(0 < f < 1) such that for some positive integer N: 

(a) d[Q{t)] c Q(f ) for < t < f < 1. 

(b) The complement C \ Cl{t) consists of N connectivity components for all t. 

(c) The boundaries dQ{t) are all C^-smooth and vary in a C^-smooth fashion with t. 

(d) The intersection Hf ^2(0 - what might be called the skeleton - is a continuum which is 
either a point or the union of a finite number of C^-smooth curves, and, in particular, 
has zero area. 

To simplify the notation, we write Gf and Pt for the Green function and the Poisson kernel of 
Cl{t), and Gf and Pf for the associated operators, respectively. 

In the setting of (a)-(d), the classical variational formula of Hadamard |5| (cf . also |9[ ITOll ) 
tells us how the Green function Gt varies with t. We prefer to state the expression in integral 
form: 



where g{Q is the rate at which the boundary ^Q(t) moves at C, e 5Q(t) along the direction of 
the exterior imit normal vector as t grows. Note that the index t is uniquely determined by 
the condition Q e (90(t). Also recall that the Poisson kernel is extended to equal outside the 
closure of the domain where it was originally defined. 

For our purposes it is important to note that, for each t, the boundaries 5Q(t) (0 < t < i) 
essentially foliate the entire domain 0(f) in the sense that their union (which is disjoint by 
assumption (a)) equals Q(f) except for the skeleton. This makes it possible to compute area 
integrals over Q(i) using the expression 



reminiscent of polar coordinates. It follows from ||3.2|| that g > almost everywhere on Q(l). 

Remark 3.1. An explicit description of a setup of the above kind and a proof of Hadamard's 
formula in that context can be found in Schippers flOll . He assumes that the boundaries dCl{t) 
are furnished by a collection of injective and non-overlapping C^-smooth homotopies Fj{t, 6) 
a = 1, . . . , N) where 6 e [0, In] (with endpoints identified). Then, at each point C = Fj{t, 6), 
a routine computation shows that g{Q = \detDFj{t,6)\/\dFj{t,d)/d6\ with DFj denoting the 
Jacobian matrix. The expression | |3.2| | now follows by usual change of variables using the maps 

3.2. The Hadamard operator. We keep the setting of (a)-(d) above, and for each C e ^2(1), 
except the points of the skeleton, let t(C) be the unique value of t such that C e dCl{t). In words, 
t(C) is the time when the front dCl{t) passes the point C- For each < f < 1, we introduce the 
Hadamard operator Qt by 



whenever the integral makes sense. In terms of signals, Qf maps the unit point mass at C £ 0.(t) 
to the Poisson kernel function z i-^ Pt(C)(z, Q, treated as zero for z e C \ Q(t(C)). Note that Qt/ 



(3.1) 




(3.2) 
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is always supported in cl[Q(f)] and Qtf = Qi(lc2(f)/) where lQ(f) stands for the characteristic 
function of Q{t). Similarly, we define the adjoint Hadamard operator 

q;/(q-1o(o(Q r p.(c)(2,Q/(z)dA(z), 

Jn(t) 

whenever the integral makes sense. 

We note the linear appearance of the Poisson kernel in the definition of the Hadamard 
operator, while it appears in a bilinear fashion in Hadamard's formula (13. 11 1. These two objects 
are intimately related; in fact, as we will see in Section|5l Hadamard's formula corresponds to 
the calculation of correlations for the stochastic field generated by the Hadamard operator. 

We collect the main properties of Qf and Q* in the following proposition. We denote by A( 
the Laplacian on Q(t), and recall that the powers of -At are to be understood with respect to 
vanishing Dirichlet boundary data. 

Proposition 3.2. The operators [-Af]-^^^Qf and QJ[-At]-^^^ are both unitary on L^(Q(f)). Moreover, 
QtQf' - [-AfAf']~^/or every t and t' , where t At' denotes the minimum oft and t'. 

In view of identity ||2.1| |, the proposition can be rephrased as saying that Qj is an isometric 
isomorphism L^(Q(t)) HQ{Q.{t)) and Qj is an isometric isomorphism H"^(Q(f)) — > L^{D.{t)), 
where H~^{D.{t)) is the dual of HQ(0(f)) under the duality pairing induced by the inner 
product. 

The following observation, which will be needed in the proof of Proposition |3.21 is helpful 
in understanding the evolution of the operator Qf as the domain Q{t) varies. 

Lemma 3.3. For < t < t' < 1 and f e L2(Q(f')), the function 

Q,/(z)-Q,/(z)= r P,(o(z,Q/(QdA(Q 

J£2(f')\£2(f) 

equals Qt>f(z) on 0.(t') \ 0.(t), and on D.(t) it is harmonic and its boundary values on dQ{t) are those 
ofQt'f. 

Proof. This follows from the easily verified fact that Qtf vanishes in C \ 0.(1), together with the 
observation that z i-^ PT(g(z, Q is harmonic in Q(t(C)). □ 

Proof of Proposition 13.21 We start by addressing the last assertion. A computation with a test 
function / e Cf^(C) gives 



QfQ;,/(z)= r (f lnit'){QPr(c){z,QP.(Q{w,Q)dA{Q]f{w)dA{w) 

= r {f r P.{z,QPr{w,Qg{Qds{QdT]f{w)dA{w), 



using Fubini's theorem and the "polar coordinates" trick i3.2\ with t as the time variable. 
However, by Hadamard's formula (13.11 1, the integral in braces equals GtAt' (z, iv). Consequently, 
QfQf'/(z) = GfAf'/(z), which gives QfQJ, = GfAc; this is equivalent to the last assertion. 
Turning to the operator Qj[-Af]^^^, we note that for all f,g & C^(Q(f)) we get 

{Q;[-At]'"f,Q;[-At]"'g)a(t) = {[-At]'"QtQ;[-At]''^f,g)Lm = {f,g)LXt) 

by what we have already obtained. Hence QJ[-Af]^^^ defines an isometry from L^(Q(f)) into 
itself, and its adjoint [-Afj^^^Qt is a contraction with dense range in L^(Q(i)). To show that both 
operators are unitary, it suffices to verify that [-Afj^^^Qf is injective. 
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So suppose that [-AtY^^Qtf = for some / e L'^{n{t)). Since Qtf belongs to Hi(Q(0), this 
implies that Qtf - 0. Then if < z < t and we replace Qtf with the function having the 
same values on 5Q(t) and being harmonic in Q(t), we still have the zero function. Thus, by 
Lemma 13. 3[ we must have Qjf = for all t with < x < f . Then the derivative must also 
vanish (for a.e. x with < x < f): 

^Qrf{z)= f P.(z,C)/(Qp(Qds(Q = 0, zeQ(x). 

But this is the harmonic extension of fg, and it can vanish everywhere only ii fg - holds 
a.e. on 5Q(x). Since g > a.e., we conclude that f - a.e. This proves that [— Af]^^^Qf is 
injective. □ 



4. Basic random fields 

We briefly recall the construction and some basic properties of the white noise and Gaussian 
free fields. We primarily adopt the viewpoint of Gaussian Hilbert spaces. See Janson's book 
Is I for an account of such spaces and Sheffield [1T| for a general survey of the Gaussian free 
field. 

Throughout this section we assume that Q is a bounded domain in the complex plane with 
a C^-smooth boundary. 

4.1. The white noise field. The (real- valued) white noise field on Q is given formally as a 
random linear combination 

i=l 

where the £,j {j = 1,2, . . .) are all independent standard Gaussian random variables, and the 
functions <pj (/ = 1,2,...) form an orthonormal basis of L^(Q). At times, we may want to express 
this as O e WN(0). The above series does not converge (almost surely) in L^(Q). However, in 
terms of the bilinear form (■, ■)□ (i.e., the inner product on Q), it makes sense to consider 

/=1 

since this series converges in mean square to a Gaussian variable with mean and variance 
||/||^2(Q) for each / e L^(Q). Thus we may view the white noise field as a Gaussian Hilbert space 
formed by the jointly Gaussian variables (/, <S>)a - one for each / e L^(Q) - whose covariance 
structure is given by 

IE[</,f)Q(^,<l')Q] = (/,^)Q. 

Here, the symbol E stands for the expectation operation. Note, in particular, that the law 
for the field is independent of the choice of basis - this is somewhat akin to the fact that 
the reproducing kernel of a Hilbert space of functions, which may be written in terms of an 
orthonormal basis, is independent of the choice of the basis (cf. [1 1). 

Alternatively, O can be imderstood as a random element in a space bigger than L^(0). One 
way is to define cD as a random distribution in the Sobolev space H~^~^{Q) for some e > such 
that its action on test functions / e C[5°(Q) yields random variables (/, <5)q as described above. 
This is possible because the natural embedding of L^(Q) into H~^~'^(Q) is a Hilbert-Schmidt 
operator; see 111. Sec. 2.2]. 
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An important property of the white noise field O is the independent action of the "vibrations" 
in different parts of Q. For instance, if Q = Oi U O2 U E, where the union is mutually disjoint, 
and Clj is open for j = 1, 2, while E has zero area, then 

1^(0) = L2(Qi)eL2(Q2) 

as orthogonal subspaces in a Hilbert space. As a result, O = Oi + O2 = ^'i ffl $2/ where Oi is 
the field O conditioned to vanish on Q2, and O2 is the field O conditioned to vanish on Qi. 
Here, the symbol ffl is used to indicate that the summands are independent random fields. In 
a natural sense, e WN(Qy) for / = 1,2. 

4.2. The Gaussian free field. The Gaussian free field on Q with vanishing boundary data, 
which we will denote as ^ e GFFo(Q), can be defined in a way that is analogous to the white 
noise field above. Instead of the space L^(Q), we just work with the Sobolev space Hj(Q) (see 
Section |Z2] |. Then we obtain a Gaussian Hilbert space formed by centred Gaussian variables 
(// ^)o,v for / e Hj(Q) with the property that 

IE[(/,W)o,v(^,^)q,v] = </,^)o,v. 

Again, W does not (almost surely) determine an element of Hq(Q), yet it can be defined as a 
random element in the Sobolev space H~^{C1) for any e > 0. In view of (|2.H , the Gaussian 
free field may be retrieved from the white noise by ^ = [— A]~^^^0 in the sense of distribution 
theory. Moreover, the action of V on L^(Q) is given by {/, W)a - ([-A]"^/, so that 

(4.1) IE[</,W)o(^,^)q] = (/,[-A]-V)o 

for f,g e L^(Q). In fact, the definition of (/, ^>o can be extended to all / e H~^{D.) since [- A]"^ 
maps H-i(Q) into Hj(Q). 

A main difference between the fields GFFo(Q) and WN(0) is that while WN(Q) is purely 
local, in GFFo(Q) we have non- trivial long-range correlations, due to the non-local nature of 
the operator [-A]"^^^. 

5. The Gaussian free field via the Hadamard operator 

Throughout this section, we work in the setting of the assumptions (a)-(d) described in 
Section |3T1 For brevity, we write Q - 0(1) for the largest domain under consideration. 

5.1. Construction. We start by invoking the Hadamard operator to produce the Gaussian free 
field from the white noise field. We agree to write Wq = 0. 

Theorem 5.1. Let O e WN(Q). For < f < 1, let Wt := QfO, i.e. 

(/,^t)o:=(Q;/,0)n, /eH-i(O). 

Then M^j e GFFo(Q(0)- Moreover, the process (0 < t < 1) has independent increments: for 
all = f(0) < ■•• < t{n) - 1 and /i, ...,/„ e H"^(Q), the random variables {fj,'^t(j) - ^f(;-i))o 
(/ = !,...,«) are independent. 

Proof. For < t,t' < 1 and /, g e H~^(Q), we have from the definitions that 

iE[(/,Wf)Q<g, Vf')o] = iE[<Q;/,o)n(Q;,g,o)c] = (q;/,q;.^)o = {f,Q,Q;g)Q. 

According to Proposition 13.21 QiQ*^, = [-AfAt']~^. Taking f = f we see that satisfies the 
correlations condition (|4.H on Cl{t). Moreover, it follows easily that the increments of are 
uncorrelated and hence, being jointly Gaussian, independent. □ 
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Since [— Af]"^ - Gt (see Sec. l2.Hl , the covariance structure of the process can be written in 
the form 

^(ma(g,'^'r)a] ^ f f Gt^r{z,zv)f{z)giw)dAiz)dA(zv) 

/c i\ Jo Jo 

VGtAt'f-'^GtM'gdA. 

In particular, the variance of (/, ^f)£2 equals the energy integral ^^IVGf/p cL4. 

In view of the definition of the Hadamard operator, we see from Theorem 15.11 that the 
Gaussian free field GFFo(Q(0) can be thought of as obtained by integrating up (in terms of the 
area integral) the harmonic fields induced by Poisson extensions of point oscillations at each 
point C e 5Q(t) with < t < t. Moreover, for < t < f < 1 the GFFq field on Q(f ) can be 
decomposed into independent fields: 

(5.2) ^Wta (W, - Wt) = QfO ffl (Q,- - Q,)0. 

Here is a GFFq on the smaller domain Q(f) with zero continuation to Q(t') \ cl[Q(f)], and 
— is a field which is harmonic in Q(i) and coincides with Wf on Cl{t') \ cl[Q{t)] (see 
Lemma l33ll . The decomposition | |5.2| | expresses the Markov property of the Gaussian free field 
(cf. El). 

5.2. Time-derivative of . Looking the increments at the infinitesimal scale, we may examine 
the time-derivative of the process ^f. To this end, it is convenient to use the harmonic sweep 
operator PJ given by 

P;/(Q:= r PKz,Q/(z)dA(z), C e 5Q(f). 
J"(f) 

It is the adjoint of the Poisson extension operator. We see from the definition of Q* that 
Q;/(C) = P;/(C) provided that C e 5Q(t) with < t < f. Consequently, for all <t,t' <1 and 
f,g & L^{(y), we have, in view of | |3.2| |, that 

^[{mn{g,^,')Q] = {Q;f,Q;,g)Li= f if p;/p;^f-ds}dT. 

Jo U50(t) > 

Writing Wt for the time-derivative of in the sense of distribution theory, we get the expression 

(5.3) ^{f,'Vt)a{g,'iJt')n]^dQ{t-t') f P;fP;ggds. 

J<?0(f) 

Next, suppose that Hf is a weighted white noise process on dD.{t) with weight g^^^, i.e., a 
Gaussian random field acting on L?-{dQ{t)) with correlation structure 

(5.4) E[((/),Sf)^,Q(,)(i/',H,)aQ(,)l = (p\}jgds 

Jda{t) 

for (/),!/'£ L^((90(f)). Let PfHf be the Poisson extension of Hf, i.e., the harmonic field on D.(t) 
given by {/, PfSf)Q(f) = (P^/, ^t}da{t) for / e L^(Q(i)); in view of | |5.4| |, it has the correlation 
structure 

(5.5) Ef</,PtH,)a(o(g,PtSt)o(J = f P;/P;^pds, 

which we compare with 1 15.31 . We arrive at the following conclusions, whose detailed verifica- 
tion is left to the reader: 
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Proposition 5.2. Suppose that (0 < t < 1) are given by Theorem \5?l\ Then the following hold: 

(a) For a fixed f , we have 

</,^f>Q = 4^)dB{T) = B{j^ K(T)dT), 0<t<l, 

where B is a standard Brownian motion with B(0) = and k(t) :- J^^-j^^^lP'/ppds. Here, "=" indicates 
that the processes have the same law (i.e., are versions of each other). 

(b) The covariance structure of the differentiated process is induced by the process PfSf.- 

^[{f,^t)a(g,^t')a] = 6o(i-i')E[</,P,S,)o(o(^,PfS,)o(o]. 
This identity should be understood in the sense of distribution theory. 

Thus, in a sense, we are dealing with a time-changed Brownian motion in the time variable t. 
On the other hand, the individual independent increments (at infinitesimal scale) correspond 
to the harmonic extension to Q(i) of weighted white noise fields along the boundaries dCl{t). 

Remark 5.3. The independent decomposition of the Gaussian free field given by Proposition |5.1| 
- and the continuous limit of that decomposition alluded to above - are in many ways analogous 
to the Brownian exploration of the Gaussian free along a space-filling curve, as in [11 1. 

5.3. Boundary averages. It is rather well known (see e.g. flT, "2]) that the averages of the 
Gaussian free field on concentric circles follow a time-changed Brownian motion. We want 
to formulate an analogue of this fact for the averages of on the associated flow of the 
boundaries dQ.{t). 

We first note that Lemma l33l implies that for / e 1^(0) and < f < 1, we have PtQif - 
(Qi - Qt)f on 0(0. Passing to the adjoints, we get Q\P*J = (Q\ - QDf for fimctions / on D.{t). 
Thus we may evaluate with the distribution P*f ds as a "test function" to obtain 

r w.rjds = (p;/,Wibow = (Q*iP;/.i')o = <(q; - q;)/,o)c = (f,w, - 

Jda(t) 

More generally, this reasoning can be extended to the case where / is replaced by a distribution 
whose support is contained in Q(i). Then the following result is obtained. 

Corollary 5.4. Let Xt{f) := J^^^^^ ^iPJ/ ds for < t < 1 and distributions f supported on the skeleton 
Pit 0(f)- Then, for all such distributions f and g, and < t,t' < 1, we have 

E[xt(/)x, ig)] = f if rj P\g Q ds] dx, 

where tvt' is the maximum of t and t'. In particular, Xt{f) = b(^J^ K(T)dT^ with jc(t) defined as in 
Propositio ns. 2\ 

Thus, the process Xt{f) {0 < t < 1) - the values of which are the integrals along the curves 
dD.(t) of the field M^i e GFFo(O) with respect to the (signed) measures PJ/ ds - obeys the law of 
a time-reversed Brownian motion running at a variable speed. 

For instance, if Cl{t) is the open disc of radius f centred at the origin and / = 6o, the unit 
point mass at 0, we have p = 1 and P^So = {2nt)~^. Then K(f) - {2nt)~^, and we obtain the 
well-known result that the average of GFFq on concentric circles of radii e~^"' has exactly the 
law of a Brownian motion. 
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Remark 5.5. It is possible to supply an expression of type | |5.1| | for the covariance ]E[Xf(/)Xt'(^)] 
altough it is slightly more delicate. We first observe that the Green potential Gi/ is well defined 
(as a distribution) for distributions / supported in the skeleton, and that G\f is harmonic in 
O off the skeleton and vanishes on dQ.. We let Gtf denote the function which equals Gif on 
Q \ Cl{t), and is harmonically extended from its boundary values to 0(f). Then 



Remark 5.6. (a) As we have seen, the GFFq field can be derived from the white noise field 
by applying the inverse of the square root of minus the Laplacian to WN. The conformal 
invariance of the Laplacian leads to the corresponding conformal invariance of GFFq. If 
we replace the Laplacian by div (which appears, e.g., in Calderon's inverse conductivity 
problem), the process retains the conformal invariance structure if we allow the weight co to 
change correspondingly. Here, the weight function co may be assumed positive and smooth. 
It should be possible to extend the use the Hadamard variation technique to this more general 
setting. We might also want to consider in the complex-valued setting the operators introduced 
by Garabedian [4J. 

(b) Hadamard's variation formula also applies in the setting of the biharmonic equation 
with vanishing Dirichlet data (see |5l for the original result, and (HQ for the integral version). 
This should lead to a corresponding Hadamard operator in this setting, and that should have 
applications to the associated random fields. 
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